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Abstract 

By using a new parametrization of the dilaton field and including a cubic term in the 
bulk scalar potential, we realize linear confinement in both meson and nucleon sectors 
within the framework of soft-wall AdS/QCD. At the same time this model also correctly 
incorporate chiral symmetry breaking. We compare our resulting mass spectra with ex- 
perimental data and find good agreement between them. 



1 Introduction 



Quantum chromodynamics has been identified as the correct theory of strong interactions 
for nearly forty years. It describes the process with large transfer momenta, like deeply 
inelastic scattering experiment, quite well. Due to the property of asymptotic freedom, 
the coupling constant is small in these cases, so the perturbation theory is reliable. As the 
characteristic energy lowers, the theory becomes strongly coupled and difficult to directly 
describe the behavior of mesons and baryons. People then develop various effective models 
exhibiting some important phenomena, such as linear confinement and chiral symmetry 
breaking, which are believed to be the features of low energy QCD. An excellent reference 
in this field is pp. 

Recently a new methodology, called "AdS/QCD" or more generally "holographic 
QCD", appears [21 E] and has been intensively studied. It is motivated by the AdS/CFT 
correspondence in string theory [U EJ |6], which identifies the M = 4 four dimensional 
supersymmetric gauge theory on one side with the type IIB string theory in AdS§ x S 5 
with Ramond-Ramond 5-form fluxes on the other side. Based on some basic setups in 
[7J E] , the authors of [9 J construct a semi- realistic holographic dual of QCD from string 
theory. By use of the general dictionary of the correspondence, AdS/QCD, the bottom- 
up counterpart of the above stringy methods, directly introduce fields corresponding to 
important operators in QCD in the bulk. This method, compared to those top-down 
constructions, has much more flexibilities to implement the phenomenological studies. 

Chiral symmetry breaking can be realized in the so-called hard-wall model where we 
only use a slice of AdS$ with an IR cutoff and the QCD scale is related to the cutoff. 
Low-lying meson spectra and couplings fit well with the experimental data. Baryons can 
also be incorporated in this model, by introducing a pair of 5d Dirac spinors for spin-| 
nucleons [10] 0, and Rarita-Schwinger fields for spin-| A resonances [13] ■ However the 
main problem in the hard-wall model is the deviation of excited state masses compared 
with data, which should grow approximately linearly due to the confinement of quarks. To 
realize linear confinement, [H] introduce a soft-wall model where the IR cutoff is removed. 
A quadratic background dilaton field is used to suppress the Kaluza-Klein mode functions 
at infinity, and gives a linear growth spectra for the vector mesons. However this model 
entangles the explicit and spontaneous breaking of chiral symmetries. This problem is 
resolved in [16] by using a modified dilaton and a quartic term in the bulk scalar potential. 
There have been lots of works on AdS/QCD. For instance, studying of meson properties 
pZ]-[25]; relation with light-front dynamics [26]; nucleons in the hard- wall model |27]-|29j. 

1 There is also a 5d skyrmion-like baryon model, see [TT1 [12] . 
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and many others. 

The aim of this paper is to realize linear confinement in both meson and nucleon sec- 
tors, and correctly incorporate chiral symmetry breaking^). We will use a new parametriza- 
tion of the background dilaton field, which is still quadratic growth at infinity to generate 
a linear growth meson spectra. While in the nucleon sector the dilaton cannot supply a 
soft-wall to discretize the spectrum, since the bulk action for Dirac spinors is first order. 
However we can obtain a linear growth nucleon spectra if the v.e.v. of the bulk scalar, 
which is the holographic dual of the chiral condensate, is asymptotically quadratic. This 
can be obtained by a cubic term in the bulk scalar potential. The remaining parts of this 
paper are organized as follows: In section 2 we will introduce the modified soft- wall model 
for meson sectors and show that how the new dilaton and the cubic term can lead to a 
linear growth meson spectra. In section 3 we describe how to realize linear confinement 
in the nucleon sector by using a new parametrization of the v.e.v of the bulk scalar. In 
section 4 we summarize our results and discuss some further issues. 



2 Meson sector 

We will consider a modified version of the soft-wall AdS/QCD model first introduced in 
[T4] . The background geometry is chosen to be 5d anti-de Sitter space 

ds 2 = G M Ndx M dx N = a 2 ( y z)(r] flu dx^dx u - dz 2 ) , < z < oo , (2.1) 

where a(z) = 1/z. Our convention of the Minkowski metric rj^ is diag(l,-l,-l,-l). There 
is also a background dilaton field &(z). The bulk action for the meson sector is 



S M = f d b xVGe-" |-^( \\F L \\ 2 + \\F R \\ 2 ) + ||M|| 2 - m 2 x \\X\\ 2 - A ||X|| 3 | 



(2.2) 



The norm of a matrix A is defined as ||>1|| 2 = Tr(A^A). Here we include a cubic term in 
the bulk scalar potential. This term allows the dilaton $(<z) and Xq(z), the v.e.v. of X, 
to be both quadratic growth as z — > oo. These are crucial for obtaining linear spectra for 
meson and nucleon sectors. We will discuss this point more carefully below. 

Other terms in ( 12.21) are as usual: Fl and Fr are the field strength of two SU{2) gauge 
fields, L and R, which are holographic duals of the left-handed and right-handed chiral 
currents. The coupling constant g\ = 12ir 2 /N c = 47r 2 is fixed by matching the two-point 
function of vector currents. The bulk scalar X, with m 2 x = —3, is dual to the chiral 



2 Thcre are already some works about this subject in different models, see e.g. [30 ) 13 l j . 
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condensate, which transforms in the bi-fundamental representation of SU(2)l x SU(2)r. 
The covariant derivative is defined as D M X = &m X — iL M X + iXR M . The trace is taken 
in the fundamental representation of SU{2). 



2.1 Bulk scalar VEV and dilaton background 

The bulk scalar X is expected to obtain a z-dependent v.e.v. as follows 

(X) = X G {z)(\ jV (2-3) 

From the bulk action (12.21) it can be seen that the possible Xq(z) and the dilaton Q(z) 
are related with each other through the EOM 

d z ( a 3 e^d z X ) - aV V^X + ^X 2 ) = . (2.4) 

The dilaton appears here because the action is second order and we need to do a integration 
by parts. This is a nonlinear differential equation for Xq, and difficult to solve for given 
dilaton profile. Here we follow the strategy of [IE], i.e. first choose a parametrization of 
Xq satisfying proper boundary conditions, then using (12. 4p to determine to background 
dilaton $ up to a integral constant 

&(z) = J^r (( « 3 ^o ') ' - a\m\X + \\Xl )) . (2.5) 

Here the prime denotes the derivative with respect to z. We assume the asymptotic 
behavior of X to be 

X Q ~ ^KCz + p), ^^0; (2.6) 

X ~ z -)• oo. (2.7) 

The normalization £ = \/3/(27r), see 02]. The parameter m 5 and <r are interpreted as 
the quark mass and chiral condensate. The second equation is needed for having a linear 
nucleon spectra, which becomes clear when we discuss the nucleon sector in the next 
section. To satisfy these two boundary conditions we choose the parametrization of X 
as follows 

„ . , Az + Bz 3 , , 

Mz) ^ vfrm ■ (2 - 8) 
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The relations between m q , a, 7 and A, B, C are 

2-4 1 ~, 2B , . 

m 9 = — , a = 2C(B--AC 2 ), 7 = • (2-9) 

By use of the relation ( 12. 5 p we can see that the asymptotic behavior of the dilaton field 
is indeed 0(z 2 ) as z — > 00. This will give us an asymptotically linear meson spectra. 
We will determine the parameters A, 5, C and A by fitting the experimental data of 
the vector and axial- vector meson masses under the constraint from the Gell-Mann-Oakes- 
Renner relation 2m q a = f 2 m 2 with m n = 139.6 Mev and f n = 92.4 Mev. The fitting 
results are: A = 0.82 Mev, B = (298.9 Mev) 3 , C = (1353.62 Mev) 2 , and 

A = 41.48. (2.10) 

Then by using (I2.9P we can obtain the values of m q , a, 7 as 

m q = 5.95 Mev, a = (241 Mev) 3 , 7 = 29.1 Mev. (2.11) 

In the next two subsections we will use them to calculate the mass spectra of vector and 
axial-vector mesons. 



2.2 Vector mesons 

From two bulk SU(2) gauge fields L M and Rm, we can introduce Vm = {L M + R M )/y/2. 
The vector p meson and its radial excitations are identified, in this model, as the Kaluza- 
Klein modes of V^. Decompose the field V^( X, Mj clS 

V^x,z) = Y,P^)f^). (2.12) 

n 

Insert this expression into the action ( 12. 2ft . with the gauge condition V z = 0, we can get 
the equation of motion for as 

- d 2 f^ + u/d M fp> = mP 2 f^ , (2.13) 

where 00 = $(2) — a' /a = $(2) + log z. After a Liouville transformation fy = e^^i/jy^, 
we can write the above equation in a Schrodinger form 

- d^P + Qu/ 2 - iu/) = mf 2 ^ . (2.14) 

The effective potential is 0(z~ 2 ) as z — > and 0(z 2 ) as z — > 00. Therefore the resulting 
spectrum is asymptotically linear. We can numerically solve this eigenvalue problem. 
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p 





1 


2 


3 


4 


5 


m ex p(Mev) 


775.5 


1465 


1720 


1909 


2149 


2265 


m t h (Mev) 


952.6 


1349 


1653 


1909 


2134 


2338 


error 


22.8% 


7.9% 


3.9% 


0.0% 


0.7% 


3.2% 



Table 1: The experimental and theoretical values for vector meson masses. The average error 
is 6.4%. 

The UV boundary condition is ipy (z — > 0) = 0. The boundary condition in deep IR 
(i.e. z — > oo) is fixed by the requirement of normalizability: /^°e~*| fv\ 2 dz < oo, with 
A being some large number. Due to the relation between fy and ipv, this is equivalent 
to | ipv\ 2 dz < oo. Since the effective potential is 0(z 2 ) at infinity, the wave function 
ipv ~ e~ Kz2 as z — > oo for some k > 0. This automatically implies d z i/jy \z — > oo) = 
used in e.g. [HI HE] . 

In Tabled] we show the predicted values of our model and corresponding experimental 
data |33j. We can see that except the ground state the agreement with the data for other 
excited states are all within 10%. The average error is 6.4%. 

2.3 Axial- vector mesons 

We define another gauge field A M = (L M — R M )/\/2. The axial- vector mesons are iden- 
tified as the Kaluza-Klein modes of A^. Decompose the field Z^j dbS 

A,( Xl z) = J2<\x)f^(z). (2.15) 

n 

Similarly choose the gauge condition A z = 0, the equation for f^\z) is 

- dlff + u/d M f£> +Agl a 2 X 2 ff = m^fP , (2.16) 

Note that, unlike the vector mesons, this axial-vector field couples to the scalar v.e.v. 
X (z), which produce a z-dependent mass term in the equation of motion. By using the 
same substitution /j^ = e^^tp^ as before, the Schrodinger equation for axial-vector 
mesons is 

- d^f + Qu/ 2 - \j< + 4 g\ a*X^j # = mf Vf> • (2-17) 

The effective potential for the axial- vector mesons is 0(z~ 2 ) as z — > and 0(z 2 ) as 
z — > oo. 
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1 


2 


3 


4 


m exp (Mev) 


1230 


1647 


1930 


2096 


2270 


m th (Mev) 


1065 


1520 


1867 


2158 


2414 


error 


13.4% 


7.7% 


3.3% 


2.9% 


6.3% 



Table 2: The experimental and theoretical values for axial-vector meson masses. The average 
error is 6.7%. 

Since the term Ag 2 a 2 XQ is also quadratic growth at infinity in our model, the asymp- 
totic slope of the axial-vector meson spectra is different from that of the vector sector. 
This fact is a direct consequence of the asymptotic behavior of the dilaton $ and the vev 
Xq, which are necessary to obtain the linear spectra in both meson and nucleon sectors. 
However it has been argued [ISl d5], based on data, that these slopes should become 
asymptotically equal. It seems to the author that, if we realize nucleon states through the 
method of [10], there would be a clash between an asymptotically linear nucleon spectra 
and asymptotically equal slopes for mesons. 

We can numerically solve this eigenvalue problem under the similar boundary condition 
ip2 \ z — >• 0) = and d z ip2 ( z ~~ °°) = 0- Since the effective potential for axial- vector 
mesons is also 0(z 2 ) at infinity, the discussion about the IR boundary condition given 
in the above subsection also applies here. In Table [2] we show the predicted values of 
our model and corresponding experimental data [33]. We can see that except the ground 
state the agreement with the data for other excited states are all within 10%. The average 
error is 6.7%. 



3 Nucleon sector 

Now we will turn to the nucleon sector. As in [10] we realize spin-^ nucleons by introducing 
a pair of 5d Dirac spinors \&i 2 which couple with the corresponding baryon operators on 
the UV boundary z = 0. The bulk action for the nucleon sector is 

S N = J d 5 xy/G e~* Tr( C K + d) , 

C, j = -g^ 1 X^f 2 - g^2X^t. (3.1) 
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Here T M = e^T A = z5%T A , and {T A ,T B } = 2r] AB with A = (a, 5). We can choose the 
representation as T A = (7°, — ry 5 ). The covariant derivatives for spinors are 

V M *1 = dM^ + ^uffEABVi-iLM*!, (3.2) 

V M *2 = d M ^ 2 + -w^IUb^ - iR M ^2 ■ (3.3) 

Here Eab = j\Pai ^b], and the nonzero components of the spin connection oj A f is u® 5 = 
— cu^ a — \ 5™. The mass m^, can be related to the dimension of the baryon operator. Its 
classical value is 5/2. Since the operator may has nonzero anomalous dimension, we can 
also treat m$ as a free parameter as suggested in [10]. This will improve our fitting result. 
The Yukawa coupling between \l/'s and X is required by the chiral symmetry breaking in 
the nucleon sector. This term will also give us a parity doublet pattern in the spectrum. 
In the nucleon sector we have two parameters, and the Yukawa coupling constant g. 
We will determine their values by fitting the nucleon mass data. 

Note that in the nucleon sector there is a crucial difference with the meson sector 
discussed in the previous section. The action for the bulk spinor is only 1st order. The 
dilaton field will not affect the equation of motion of bulk spinors, since no integration by 
parts is needed here. However we will see explicitly below that their Yukawa couplings 
with the bulk scalar will give us a quadratic growth soft- wall if the v.e.v X (z) ~ 0(z 2 ) 
as required in (12. 7p . 



3.1 Nucleon spectra 

As in the above meson sectors, nucleons in this model are identified as the Kaluza-Klein 
modes of the 5d Dirac spinors 

&i(x,z) - , , , , ; * 2 foz) , , , , • (3.4) 

Here N^ 1 ' are two-component Weyl spinors, and are complex scalar functions. The 
internal components should satisfy [10] 

-mf (J I , (3.5) 

\J2R 

+ m<?> I J I . (3.6) 

, J2L 
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Here A^ 1 = 2 ± and 

v(z) =gX {z) a{z) = • (3.7) 

V 1 + C 2 z 2 

To get the correct chiral coupling between \l/'s and baryon operators at the UV boundary, 
we need to impose boundary conditions [10] 

A ( ^->0) = 0, /S(^0) = 0. (3.8) 

For the boundary condition at infinity, we choose 

/$(*-► oo) = 0, fS(z ^oo) = (3.9) 

to guarantee the normalizability of eigenfunctions. We can eliminate (/ijs\/2b) T fr° m 
EOM ([33D and ($3M) as 




(n)2 m\ +171^—6 o / 
-l/ 



+ _ . J =0.(3.10) 



N 




Here u' means d z v. We can transform the above equation to a coupled Schrodinger form. 
Define (fg,f£>) = (z 2 ^, z\S) and X J*> = (xS,xSf, then 

Xi n) " + (m& )2 -^))xi n) = 0, (3.11) 

where the potential matrix V(z) is 





v 12 \ 


H 


\V 21 


V 22 ) 





m * _m >T, l 2 / 
Z SL _L n z 11' 



Note that the diagonal entries Vii and V22 are both 0(^ 2 ) as 2; — > 00 due to the v 2 
term, and the off-diagonal ones tend to a constant. Therefore the nucleons, like mesons 
discussed in the previous section, also have an asymptotically linear spectrum. This is 
why we require X (z) should has the behavior as in (12. 7p . The boundary condition for 
Xl } can be deduced from that of in (ESJ), (JUS} and the EOM 1ET5]) . (13T61) as follows 

xi? = 0, ^xS + xS' = 0, atz^O; (3.13) 
xg = 0, ^X$-xg' = 0, atz^oo. (3.14) 



N 





1 


2 


3 


4 


5 


6 


m exp (Mev) 


940 


1440 


1535 


1650 


1710 


2090 


2100 


m th (Mev) 


1144 


1398 


1509 


1698 


1798 


1955 


2046 


error 


21.7% 


2.9% 


1.7% 


2.9% 


5.2% 


6.5% 


2.6% 



Table 3: The experimental and theoretical values of the spin-1/2 nucleon masses in the one- 
parameter fit. The average error is 6.2%. 



N 





1 


2 


3 


4 


5 


6 


m exp (Mev) 


940 


1440 


1535 


1650 


1710 


2090 


2100 


mth(Mev) 


1003 


1333 


1477 


1701 


1826 


2004 


2117 


error 


6.7% 


7.4% 


3.8% 


3.1% 


6.8% 


4.1% 


0.8% 



Table 4: The experimental and theoretical values of the spin-1/2 nucleon masses in the two- 
parameter fit. The average error is 4.7%. 

We use the values of A, B, C or m q , a, 7 coming from the meson sector. As a first 
step, we fix m# to its classical value 5/2 which is determined by the classical dimension 
of the baryon operator, and only use the parameter g to fit seven spin-| nucleon masses 
listed in PDG [33]. The best value is 

# = 16.205. (3.15) 

The resulting mass spectra are listed in Table |3j We can see from this table that even 
this one-parameter fit is not bad. Except the ground state the agreement with the data 
for other excited states are all within 10%. The average error is only 6.2%. 

Since the corresponding baryon operator may has nonzero anomalous dimension, we 
can treat the mass also as a free parameter. By fitting the experimental data, the 
preferred values are 

m*=1.5, 3 = 19.515. (3.16) 

In Table H] we show the predicted values of our soft-wall nucleon model and the correspond- 
ing experimental data. We can see that this two-parameter fit improves the agreement. 
Even the worst one is smaller than 7.5%. The average error is only 4.7%. 

4 Summary and discussion 

In this paper we introduce a modified soft- wall AdS/QCD model in which linear confine- 
ment is realized in both meson and nucleon sectors. Linear confinement in the meson 
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sector require the dilaton to be asymptotically quadratic at infinity, while in the nucleon 
sector it requires the v.e.v. of the bulk scalar to be quadratic growth. The crucial point for 
these two conditions being satisfied at the same time is the cubic interaction term in the 
bulk scalar potential^ We also compare the predicted mass spectra with the experimental 
data and find good agreement between them. 

There are also some further issues, beyond the scope of this paper, need to be studied 
in the future. For example, extend this model to spin-| A resonances. More importantly, 
we could study the coupling between mesons and nucleons. In the hard-wall model it 
seems that mesons and nucleons need different IR cut-offs to fit the data. This problem 
becomes much more serious when we consider their couplings. However in our soft- 
wall model presented in this paper, this consistency problem disappears by definition. 
Therefore it seems more natural to study the meson-nucleon couplings in our model. 
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